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TECHNICAL  ROTE  2606 

SPECTRUM  OF  TURBULENCE  IN  A  CONTRACTING-  STREAM 
By  H.  S.  Ribner  and  M0  Tucker 


SUMMARY 

The  spectrum  concept  is  employed  to  study  the  selective  effect  of 
a  stream  contraction  on  the  longitudinal  and  lateral  turbulent  velocity 
fluctuations  of  the  stream.  By  a  consideration  of  the  effect  of  the 
stream  contraction  on  a  single  plane  sinusoidal  disturbance  wave, 
Mathematically  not  dissimilar  to  a  triply-periodic  disturbance  treated 
by  G.  I.  Taylor,  the  effect  on  the  spectrum  tensor  of  the  turbulence 
and  hence  on  the  correlation  tensor  are  determined.  Lack  of  inter¬ 
ference  between  waves  follows  from  the  postulation  of  a  low  level  of 
turbulence;  this  and  the  assumption  of  an  inviscid  fluid  imply  neglect 
of  decay  effects.  The  compressibility  of  the  main  stream  is  taken 
into  account,  but  the  density  fluctuations  associated  with  the  tur¬ 
bulence  is  assumed  to  be  negligible;  this  would  be  the  case  if  the 
turbulence  originated  from  wakes  and  boundary  layers  in  the  very  low 
speed  portion  of  the  flow.  For  an  axisymmetric  contraction  and  a 
particular  isotropic  initial  turbulence  some  explicit  results  are 
obtained.  The  one -dimensional  longitudinal  spectrum  is  found  to  be 
distorted  (as  well  as  reduced  in  amplitude)  with  its  peak  shifted  well 
to  the  right  of  the  initial  position  above  the  zero  of  the  wave-number 
scale.  The  selective  effect  of  the  contraction  on  the  mean  square 
longitudinal  and  lateral  components  of  turbulent  velocity  is  found  to 
be  given  uniquely  when  the  initial  turbulence  is  isotropic,  regardless 
of  the  details  of  the  spectrum.  If  the  initial  spectrum  is  anisotropic 
as,  for  instance,  that  produced  by  a  damping  screen,  then  the  selective 
effect  is  altered. 

In  a  crude  extension,  decay  effects  outside  the  scope  of  the 
theory  are  allowed  for  in  first  approximation.  With  this  extension, 
a  comparison  with  experiment  is  made  of  the  selective  effect  on  tur¬ 
bulent  intensity  where  the  estimated  decay  effects  are  comparable  with 
the  contraction  effects.  The  agreement  is  good  for  the  longitudinal 
component,  very  poor  for  the  lateral  component,  the  experimental  data 
themselves  being  in  conflict. 
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OTTEODUCTION 

The  generation  of  a  wind-tunnel  flow  is  always  accompanied  by  a 
certain  amount  of  turbulence;  this  is  one  respect  in  which  the  flow  fails 
to  simulate  free-f light  conditions.  Measurements  in  the  tunnel,  particu¬ 
larly  those  sensitive  to  boundary-layer  behavior,  are  known  to  be  affected 
by  this  turbulence.  Accordingly,  the  tunnel  designer  attempts  to  reduce 
the  intensity  to  the  lowest  practicable  level.  The  use  of  honeycombs  and 
damping  screens  in  a  large  low-speed  section  (settling  chamber)  followed 
by  a  sharp  contraction  to  the  much-higher-speed  working  section  is  known 
to  be  effective.  The  honeycombs  and  screens  located  in  a  low-speed 
section  reduce  the  absolute  level  of  the  turbulence  with  little  drag 
penalty;  then  the  relative  level  is  greatly  reduced  by  the  large  gain 
in  tunnel  speed  through  the  contraction,  aside  from  any  effect  of  the 
contraction  on  the  absolute  level. 

Once  the  characteristics  of  honeycombs  and  screens  are  known,  the 
further  quantitative  estimate  of  the  reduction  in  turbulence  involves  a 
knowledge  of  the  effect  of  the  tunnel  contraction1  on  the  turbulence. 

It  is  known  that  the  longitudinal  component  of  the  turbulence  is  greatly 
reduced  (in  absolute  value)  by  the  contraction;  the  behavior  of  the 
lateral  component  appears,  on  the  other  hand,  to  vary  from  no  change  to 
a  substantial  increase.  Prandtl  (reference  l)  obtained  a  quantitative 
estimate  of  the  first  effect  by  considering  the  conservation  of  energy 
for  a  perturbed  longitudinal  filament:  if  the  initial  stream  speed 
is  U,  the  filament  speed  U  +  u,  with  u<<U,  and  the  final  stream  speed 
is  7U,  then  the  final  filament  speed  must  be  2U  +  2“^u;  that  is,  the 
contraction  reduces  the  longitudinal  perturbation  velocity  u  by  the 
factor  Z"*1.  For  the  lateral  effect,  Prandtl  applied  conservation  of 
momentum  to  a  small  rotating  cylinder  of  the  fluid,  with  its  axis 
cross  stream,  as  the  fluid  traversed  the  tunnel  contraction.  He  con¬ 
cluded  that  the  lateral  perturbation  velocity  v  is  increased  by  a 
factor  -y/T. 

Prandtl* s  considerations  on  the  effect  of  a  stream  contraction 
were  limited,  as  has  been  noted,  to  particular  idealized  "eddies". 

G-.  I.  Taylor  later  (reference  2)  attempted  more  realism  by  treating  a 
mathematically  defined  model  of  turbulence  which  amounted  to  vortices 
in  parallelepiped  partitions  arranged  in  a  regular  three-dimensional 


4he  considerations  of  this  paper  are  not  limited  to  a  wind-tunnel 
contraction:  they  may  be  applied  to  any  stream  tube  of  varying  cross 
section  large  compared  with  the  scale  of  the  turbulence. 
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array.  The  changes  In  vorticity  on  traversing  the  contraction  were 
determined  from,  a  theorem  based  on  conservation  of  circulation  for  an 
inviscid  fluid;  the  corresponding  altered  turbulent  velocity  pattern 
was  then  calculated.  The  final  result  of  the  analysis  consisted  in 
expressions  for  the  root  mean  square  longitudinal  and  lateral  turbulent 
velocity  components  u’  and  v'  downstream  of  the  contraction  expressed 
as  ratios  of  the  corresponding  values  upstream. 

The  initial  condition  of  isotropic  turbulence  (mean  values  unaffected 
by  rotation  or  reflection  of  axes)  was  approximated  by  specifying  the 
vortex  partitions  to  be  cubical.  For  this  case  the  reduction  in  the 
longitudinal  component  u’  was  found  to  vary  more  nearly  like  1.5  Z-1 
than  the  value  2"1  suggested  by  Prandtl.  Wo  explicit  result  was 
found  for  the  variation  of  the  lateral  component,  however:  the  calcu¬ 
lations  contained  a  free  parameter. 

Taylor's  results  for  the  longitudinal  component  agreed  fairly  well 
with  the  experimental  data  then  available,  but  it  is  now  considered  that 
the  measurements  were  made  too  close  behind  the  screens  for  the  screen- 
produced  turbulence  to  have  been  isotropic.  On  theoretical  grounds,  the 
objection  to  Taylor’s  theory  is  threefold:  first,  the  decay  processes 
of  turbulent  mixing  and  viscous  dissipation,  which  result  in  a  reduction 
of  the  mean  intensity  with  axial  distance  in  the  wind  tunnel,  are  neg¬ 
lected;  second,  the  assumed  model  of  turbulence  fails  to  exhibit  the 
spatial  and  temporal  randomness  of  actual  turbulence;  third,  no  choice 
of  the  parameters  in  Taylor's  model  corresponds  to  isotropy.  In  a  sense 
all  three  objections  apply  likewise  to  Prandtl’ s  results:  no  model  was 
employed  in  his  considerations,  and  hence  no  distinctions  between  the 
effects  of  isotropy  and  anisotropy  were  made. 

The  second  and  third  objections  can  be  removed  by  working,  not  with 
a  model  of  turbulence,  but  instead  with  a  Fourier  integral  representation 
of  a  random  turbulent  field.  The  integral  represents  a  superposition  of 
plane  transverse  sinusoidal  waves  of  all  wave  lengths  and  with  apparently 
random  phases  and  planes  of  polarization.  This  aggregate  of  plane  waves 
constitutes  the  (three-dimensional)  spectrum  of  the  turbulence.  Only  the 
statistical  aspects  of  this  spectrum  will  be  known,  not,  for  example, 
the  detailed  phase  relationships.  Mean  square  velocity  components  may 
be  obtained  by  an  integral  of  certain  spectrum  functions  in  which  the 
phase  relationships  are  suppressed;  these  functions  are  included  in 
the  "spectrum  tensor"  (reference  3). 

Taylor’s  concepts  may  be  applied  to  find  the  effect  of  a  stream 
contraction  on  a  single  plane  wave.  The  effect  under  the  assumptions 
is  linear;  therefore  the  superposition  implied  by  the  Fourier  integral 
may  be  employed  to  obtain  the  contraction  effect  on  a  field  of  turbu¬ 
lence.  In  particular,  if  the  initial  spectrum  tensor  is  known  the  final 
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spectrum  tensor  is  determined;  the  initial  and  final  mean  square  velocity 
components  then  result  from  quadratures.  Indeed,  from  the  same  informa¬ 
tion,  changes  due  to  the  stream  contraction  in  correlations  of  velocity 
at  different  points  may  he  calculated:  use  is  made  of  the  fact  that  the 
correlation  tensor  is  an  inverse  Fourier  transform  of  the  spectrum 
tensor  (reference  3). 

Accordingly,  the  injection  of  the  spectrum  point  of  view  into 
Taylor fs  original  concept  of  the  contraction  effect  makes  possible  a 
more  realistic  calculation  of  the  changes  in  mean  square  velocity  com¬ 
ponents  .  In  addition,  it  provides  much  more  detailed  information  con¬ 
cerning  changes  in  the  statistical  properties  of  the  turbulence;  that 
is,  it  provides  the  changes  in  the  spectrum  tensor  and  in  the  correla¬ 
tion  tensor. 

The  ideas  just  outlined  are  developed  in  the  present  paper.  The 
first  section  is  devoted  to  an  account  of  turbulent  spectrum  analysis 
in  a  form  specially  adapted  to  the  analysis  of  the  contraction  effect c 
In  this  account,  which  is  a  generalization  of  a  development  in  refer¬ 
ence  4,  the  role  of  the  spectrum  tensor  is  subordinated  to  that  of  the 
individual  Fourier  components  (plane  waves)  in  contradistinction  to 
the  customary  treatment.  This  approach  has  perhaps  an  auxiliary  merit 
in  providing  some  better  physical  insight  into  the  significance  of  the 
spectrum  tensor. 

Next  the  effect  of  a  stream  contraction  on  a  single  plane  wave  is 
calculated  by  an  application  of  Taylor* s  concepts.  The  treatment  is 
slightly  more  general  in  that  compressibility  of  the  main  stream  is 
allowed  for.  The  density  fluctuations  associated  with  the  turbulence 
are  assumed  to  be  negligible;  this  would  be  the  case  if  the  turbulence 
originated  entirely  from  boundary  layers  and  wakes  in  the  very  low 
speed  portion  of  the  flow.  Following  Taylor,  the  problem  is  linearized 
by  postulating  a  sufficiently  weak  turbulence  so  that  the  self -distortion 
of  the  turbulent  eddies  is  small  compared  with  the  distortion  Imposed 
by  the  contraction  of  the  main  stream;  this  together  with  the  assumption 
of  an  inviscid  fluid  implies  neglect  of  the  decay  of  the  turbulence. 

In  succeeding  sections  the  spectrum  and  correlation  tensors  down¬ 
stream  of  the  contraction  are  expressed  in  terms  of  the  corresponding 
initial  tensors.  For  the  special  case  of  an  axisymmetric  contraction 
and  isotropic  initial  turbulence  the  ratios  of  the  root  mean  square 
longitudinal  and  lateral  velocity  fluctuations  downstream  and  upstream 
are  obtained  explicitly  in  terms  of  the  parameters  defining  the  con¬ 
traction.  For  a  particular  subcase  where  the  initial  isotropic  spectrum 
tensor  is  specified,  the  corresponding  *  one -dimensional 1  spectrums  .(as 
would  be  recorded  by  stationary  hot-wire  probes)  upstream  and  downstream 
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of  the  contraction  are  calculated;  the  specification  is  such  that  the 
upstream  one -dimensional  spectrum  corresponds  to  experiment  (refer¬ 
ence  5)  in  a  number  of  cases  of  isotropic  turbulence. 

Most  of  the  calculated  contraction  effects  are  amenable  to  experi¬ 
mental  checks  either  directly  or  indirectly.  The  available  experimental 
data,  however,  are  limited  to  the  changes  in  the  root  mean  square  veloc¬ 
ities.  A  comparison  with  these  experimental  data  is  given  with  an  esti¬ 
mated  allowance  for  decay  effects  outside  the  scope  of  the  theory.  Design 
curves  of  the  changes  in  the  root  mean  square  velocity  components  neg¬ 
lecting  decay  are  included  for  engineering  purposes. 


SPECTRUM  ANALYSIS 

Representation  of  turbulence  by  superposition  of  plane  sinusoidal 
waves.  -  Suppose  q.j_,  q.g,  q3  represent  the  components  of  velocity  in  a 

turbulent  field;  that  is,  q.^,  q.g,  and  q^  vary  in  an  apparently  random 
manner  in  space  and  time,  and  the  mean  values  q.p  ^  s  ^3  ~  0-  Sub¬ 
ject  to  certain  conditions,  a  snapshot  of  this  field  at  any  instant  can 
be  represented  as  a  set  of  three-dimensional  Fourier  integrals 


<la(Xi,  Z2,  X3)  = 


i(k1x1+k2x2+]£3x3) 


dk1dk2dk3; 


where  a  =  1,  2,  or  3  and  the  significance  of  k-^,  kg,  and  k^  will  be 

brought  out  later.  A  continuous  representation  of  the  turbulent  field 
is  obtained  by  allowing  the  Qa  to  vary  with  time. 


It  will  be  convenient  to  abbreviate  the  Fourier  integral  to 


(la) 


and  to  introduce  the  companion  equation 


l 
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where  k  =  k]_,  k2,  k3;  x  =  x-^,  x2,  x3;  dx (x)  =  dx-]dx2dx3.  The  second 

equation  allows ,  in  principle  at  least,  the  coefficients  Qa(k)  to  he 

calculated.  Mathematically,  qa(x)  and  Qa(k)  are  termed  three- 

dimensional  Fourier  transforms  of  each  other;  use  will  he  made  of  this 
relationship  later. 

The  velocity  components  qa  are  connected  hy  the  condition  of  con¬ 
tinuity.  In  many  cases  of  practical  interest  these  turbulent  velocities 
originate  from  boundary  layers  and  the  wakes  of  obstacles  in  flows  of 
low  subsonic  speed,  so  that  associated  density  fluctuations  may  be 
ignored;  this  is  still  permissible  when  the  turbulence  so  produced  is 
transported  by  a  high-speed  stream.  Thus  the  incompressible  form  of 
the  continuity  equation  may  be  used  and  the  result  is 

Q^_k]_  +  Q2k2  +  Q3k3  #=  0 
This  relation  may  be  written  more  compactly  as 

=  0  (2) 

a 

Physical  interpretation.  -  The  amplitude  components  Qa  are  com¬ 
plex  in  general.  According  to  equation  (lb),  then,  the  requirement  that 
the  velocity  components  qa  be  real  implies  that  Qa(-k)  is  the  complex 

conjugate  of  Q^k) .  If  corresponding  terms  for  k  and  -k  in  equa¬ 
tion  (la)  are  paired  their  sum  is  thus  equal  to  the  real  quantity 


2 (Re  Qa)  cos  (k'x)  -  2(lm  Qa)  sin  (k*x) 


(3) 


The  imaginary  parts  cancel  in  the  pairing,  which  implies  that  they  con¬ 
tribute  nothing  to  the  integral.  Expression  (3)  represents  a  pair  of 
plane  standing  waves,  a  cosine  wave  and  a  sine  wave,  with  normals  in 
the  direction  k  =  (k-j_,  k2,  k3),  where  x  =  (x1;  x2,  x3)  is  the  radius 

vector  to  any  point.  The  vector  k  is  termed  the  wave-number  vector 
and  its  magnitude  k  simply  the  wave  number;  the  corresponding  wave 
length  is  2n  divided  by  the  wave  number.  Since  k  is  perpendicular 
to  the  wave  front,  it  is  sometimes  referred  to  herein  as  the  'wave 
normal ' . 

The  continuity  condition,  equation  (2),  states  that  both  the  real 
part  (Re  Qa)  and  the  imaginary  part  (im  Qa)  of  the  amplitude  vector 
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Q  =  Qa  =  (Q^,  Q2;  Qg)  are  perpendicular  to  the  wave  normal  k;  that  is, 

both  -waves  of  expression  (3)  are  transverse.  For  each  wave  any*  one  of 
the  parallel  planes  containing  both  the  local  velocity  vector  q  and 
the  wave -number  vector  k  (which  is  perpendicular  to  Q  and  hence  to 
q)  is  called  the  plane  of  polarization0  The  cosine  wave  (real  part) 
and  sine  wave  (imaginary  part)  may  be  polarized  in  different  planes  in 
general;  the  necessary  and  sufficient  condition  they  be  polarized  in 
the  same  plane  is  that 


Re  Re  Q2  Re  Q3 

Im  ~  Im  Q2  ”  Im  Qg  ^ 

Equations  (lb)  are  now  seen  to  represent  a  superposition  of  plane 
sinusoidal  waves  (Fourier  components)  with  all  orientations  of  the  wave- 
fronts  (all  directions  of  the  wave-normal  k)  and  all  wave  lengths  (all 
wave  numbers  k).  Each  wave  is  transverse,  and  all  planes  of  polariza¬ 
tion  are  permitted 0  For  each  value  of  k  there  exists  a  cosine  wave 
and  a  sine  wave;  their  respective  amplitudes  and  planes  of  polarization 
are  different  in  general.  The  complex  amplitude  components  Qa(k) 

express,  in  their  real  and  imaginary  parts,  how  the  respective  ampli¬ 
tudes  and  planes  of  polarization  vary  with  the  wave-front  orientation 
and  the  wave  number. 

Mean  values:  the  correlation  tensor.  -  Consider  the  spatial"*" 
mean  value  of  the  product  of  the  velocity  component  qa  at  x  and  the 

velocity  component  q^1  at  x1  =  x  +  r  as  x  varies  but  the  separation 

r  of  the  two  points  remains  fixed  during  the  averaging  process;  this 
mean  value  is  called  a  velocity  correlation  and  is  given  the  symbol 
Rap(r).  There  are  nine  such  correlations,  corresponding  to  a  *  1,  2,  3, 

P  =  1,  2,  3.  The  form  R^p(r)  has  been  shown  to  transform  like  a 

second-order  tensor  and  has  been  designated  (sometimes  divided  by  q^) 
as  the  "correlation  tensor"  (reference  6). 


-*-If  the  statistical  properties  of  the  turbulence  are  independent  of 
position  (homogeneous  turbulence)  and  time-independent,  an  average  at  a 
given  time  over  all  space  equals  an  average  at  a  given  point  (or  pair  of 
correlated  points)  over  all  time;  a  proof  is  given  in  appendix  B.  If  the 
statistical  properties  vary  slowly  with  time  the  space  average  will  still 
approximate  a  time  average  over  an  interval  just  long  enough  to  smooth 
out  the  fluctuations. 
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form  rap(k)  defined  by- Batchelor  (reference  3)  aa  the  Fourier  transform 
of  Rap(r) .  The  form  ragk)  is  known  as  the  spectrum  tensor.  The 

Fourier  transform  relations  connecting  the  spectrum  tensor  and  the  corre¬ 
lation  tensor  are  summarized  as 


(6a) 


*043(21) 


(6b) 


By  use  of  the  Fourier  transform  relationship ,  Batchelor  demonstrated 
that  rap  is  a  second -order  tensor  and  obtained  a  number  of  its  proper¬ 
ties.  Thus,  for  example,  Tap  is  complex,  in  general,  with  Tpa  =  1^*, 
and  the  diagonal  elements  rnrr  are  real;  also,  fap(-k)  m  rpa(k) . 


It  is  of  interest  to  observe  that  these  same  properties  result  immedi- 

XjXul  0 

ately  from  the  identification  of  T_>.0O  QaQp*  with  rap.  Thus 
QaQp*  is  complex,  in  general;  QpQo,*  equals  [QaQp*]*^  ^  ia> 

of  course,  real.  Furthermore,  since  Qg-k)  =  Qa*(k),  Qa(-k)Qp*(-k) 
equals  Qc*(k)Qp(k);  hence  ^(-k)  =  rpa(k). 

The  distinction  "between  cases  where  is  real  and  cases  where 

it  is  complex  may  "be  given  a  physical  interpretation.  The  product 
QaQp*,  and  correspondingly  rap,  is  seen  to  "be  real  when  the  condition 

equation  (4)  is  satisfied.  This  implies  that  the  cosine  wave  and  sine 
wave  associated  with  wave  number  k  are  polarized  in  the  same  plane. 

The  alternative  condition  where  and  hence  are  complex 

implies  polarization  of  corresponding  cosine  and  sine  waves  in  different 
planes.  The  velocity  pattern  of  such  a  pair  of  waves  is  quite  inter¬ 
esting:  successive  velocity  vectors  along  a  line  in  the  direction  of 
the  wave  normal  k  turn  progressively  about  this  line  in  spiral  fashion; 
the  tips  of  the  vectors  trace  out  a  helical  curve  on  a  cylinder  of  oval 
cross  section. 


Energy  spectral  density.  -  Each  of  the  diagonal  elements  T -q, 
V22,  and  r33  of  the  spectrum  tensor  Fap  may  "be  interpreted  as  an 
energy  spectral  density.  Thus,  according  to  equation  (6b) 


« 


V 


* 


p 


4 
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%l(0)  s  ux  = 


rn(k)dT(k);  therefore 


the  differential 

energy  component 
within  the  range 


|rudt  (k)  represents  the  contribution  to  the  kinetic 

—  U]_2  per  unit  mass  made  by  waves  with  wave  number 
2 

dT(k)  o 


One -dimensional  spectrum.  -  The  elements  of  the  three-dimensional 
spectrum  tensor  are  not  directly  measurable;  they  may  he  obtained  by 
taking  the  Fourier  transform  of  the  measured  correlation  tensor.  A  hot¬ 
wire  probe  placed  in  the  moving  stream  will,  however,  develop  a  fluctu¬ 
ating  output  voltage  whose  (one -dimensional;  frequency  spectrum  (refer¬ 
ence  7)  is  related  to  a  diagonal  element  of  the  three-dimensional  spec¬ 
trum  tensor.  Thus  by  equation  (6b)  the  contribution  to  the  mean  square 

velocity  component  q  2  (=Earr  (0) )  from  all  waves  with  wave-number  com¬ 
ponents  in  the  direction  of  the  x-^axis  between  j  k-^j  and  jk^j  +  JdkJ 
is 

Fa(ki)dki=  2 

the  factor  of  2  accounting  for  suppression  of .negative  values  of  k^. 

The  function  Fa(k-^)  is  the  one -dimensional  spectrum  corresponding  to 
the  velocity  component  qa;  the  values  a  =  1,  2,  3  correspond  respec¬ 
tively  to  the  longitudinal  and  two  lateral  spectrums.  The  particular 
spectrum  obtained  depends  on  the  arrangement  of  the  hot-wire  probe 
elements. 


(7) 


EFFECT  OF  STREAM  COMPACTION 

Consider  now  that  the  turbulent  velocity  pattern  q^,  q^,  q^  is 
carried  along  by  an  inviscid  general  stream  with  velocity  TJ(x^)  in  the 
x^-direction.  Consider  also  that  q^,  q2,  q^  are  so  small  that  their 

effect  on  the  streamlines  may  be  neglected  as  the  flow  traverses  a  wind- 
tunnel  contraction.  The  contraction  will,  however,  distort  the  shape 
of  fluid  elements.  (See  fig,  1.)  The  vorticity  distribution  will  be 
forced  to  alter  accordingly  to  conserve  the  circulation  about  each 
element.  The  net  result  will  be  an  altered  pattern  of  turbulence.  Each 
plane  wave  (Fourier  component)  Q-^,  Q,2j  Q^e  ~  2:  will,  in  fact,  be 


12 


MCA  TN  2606 


altered  independently  under  the  linearizing  assumption  to  be  made;  the 
over-all  effect  on  q2,  q.3  can  be  obtained  by  the  summation  expressed 

by  equation  (l).  It  thus  suffices  to  consider  the  effect  of  the  contrac¬ 
tion  on  a  single  representative  plane  wave* 


Effect  of  Contraction  on  Eepresentative  Plane  Wave 

Velocity  and  vorticity  at  upstream  station.  -  Designate  by  A  a 
reference  station  upstream  of  the  contraction  and  by  B  a  reference 
station  downstream  of  the  contraction,  (See  fig.  1.)  Let  a  typical 
Fourier  component  (plane  wave)  of  the  turbulent  field  qa  (a  »  1,  2,  3) 

at  station  A  be  represented  at  time  t  =  0  by 

3/  '  Ce1<£'s)  <8) 

This  wave,  equation  (8),  is  supposed  to  he  carried  along  by  the  main 
stream,  with  velocity  U0 

N 

The  vorticity  is  obtained  from  the  curl  of  equation  (8)  as 


E‘aSrVA1(£-’ 

P>7 


where 


if  any  pair  of  subscripts  are  equal 
if  0,(3/  are  in  cyclic  order 
if  0(37  are  in  anticyclic  order 


O) 


Distortion  of  a  fluid  element  in  passage  through  contraction.  - 
Suppose  the  contraction  is  such  that  the  stream  velocity  U  is  increased 
hy  a  factor  Z^  between  stations  A  and  B  while  the  breadth  and 

height  of  the  tunnel  are  reduced  by  factors  Z2  an<^  respectively. 

(See  fig.  1(a).)  In  traveling  from  A  to  B  an  initially  cubical 
element  of  fluid  of  edge  D  will  be  distorted  into  a  parallelepiped  of 
edges  Z]D,  Z2D,  Z3D  (see  fig.  1(b));  a  particle  in  the  element  originally 

(t  =  0)  a  vector  distance  x  from  a  comer  particle  will  finally  (t  =  t) 
be  found  a  distance  £  from  the  comer  particle,  where  £  is  related 
to  x  by 
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=  l1x1 

S  2  -  ^2x2 
5 3  =  ^3x3 


(10) 


In  this  argument  the  modification  of  the  streamlines  due  to  the  turbulent 
velocity  fluctuations  has  been  neglected.  This  implies  that  the  relative 
displacement  of  two  adjacent  particles  due  to  the  superposed  turbulent 
motion  is  small  compared  with  the  displacement  due  to  the  tunnel  contrac¬ 
tion.  This  key  assumption,  due  to  Taylor  (reference  2),  linearizes  and 
vastly  simplifies  the  problem.  The  limitations  Imposed  by  the  assumption 
are  discussed  later  under  "Decay  Considerations”. 

The  velocity  ratio  2^  and  the  lateral  and  vertical  contraction 
ratios  l ^  and  l ^  are  related  by  the  continuity  condition 


°^1^2^3  ~  ^ 

where  a  is  the  ratio  of  stream  densities  at  stations  B  and  A;  the 
density  is  considered  uniform  at  each  station  in  accordance  with  the 
initial  assumption  of  negligible  turbulent  density  fluctuations. 

Yorticity  at  downstream  station.  -  The  vorticity  is  carried  along 
by  the  flow,  the  fluid  elements  undergoing  the  distortion  pictured  in 
figure  1(b),  to  the  approximation  used.  During  the  motion  the  strength 
changes  in  such  a  way  as  to  maintain  the  constancy  of  circulation  of  the 
fluid  elements.  The  changes  are  expressed  by  the  equations  for  the 

transport  of  vorticity  in  the  Lagrangian  form,  due  to  Cauchy  (see 
reference  8), 


^This  statement  is  exact  for  the  postulated  inviscid  fluid .  The 
modification  produced  by  the  diffusive  effect  of  viscosity,  in  the  case 
of  a  gas,  becomes  appreciable  for  the  smaller  eddies  or  higher  wave 
numbers;  for  this  analysis  a  criterion  for  neglect  of  viscous  effects  is 

vk2<< |du/dx|,  (See  "Decay  Considerations",  equation  (43).) 
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where  a  is  the  density  ratio  between  stations  B  and  A,  and  the  deriv¬ 
atives  c)€  /SxR  express  the  effect  of  the  fluid  distortion.1  Evaluation 

CO  p 

by  means  of  the  distortion  equations  (10)  yields  simply 


or,  in  expanded  form 


(11a) 


These  equations  relating  downstream  and  upstream  vorticity  embody  the 
entire  dynamics  of  the  contraction  effect.  The  equations  are  not  limited 
to  the  plane  sinusoidal  waves  discussed  earlier,  hut  apply  to  any  (weak) 
vorticity  distribution  whatsoever. 

The  above  derivation  of  the  vorticity  changes  is  substantially  in 
the  form  given  originally  by  G.  I.  Taylor  (reference  2)  for  the  case 
o-l  (incompressible  flow).  In  order  to  assess  the  influence  of  the 
simplifying  assumptions  a  more  general  derivation  based  on  the  Navier- 
Stokes  equations  is  given  in  the  section  entitled  "Decay  Considerations". 

By  virtue  of  equation  (9)  as  applied  to  (lla)  the  vorticity  at 
station  B  is  obtained  explicitly  as 


4 

r 


eap7^:p^'ye 


ik*x 


(11b) 


where,  it  will  be  remembered,  x  is  the  radius  vector  to  a  fluid  particle 
at  time  t  =  0  when  the  fluid  element  is  at  station  A  in  the  moving 
coordinate  system  of  figure  1(b).  The  corresponding  vector  to  the 
particle  at  time  t  =  t,  when  the  fluid  element  is  at  station  B,  is  £ 
in  that  figure.  When  equations  (10)  are  used  to  express  x  in  terms 
of  £,  the  exponential  k-x  becomes,  in  expanded  form, 


^These  equations  refer  to  axes  moving  with  some  fluid  particle 
rather  than  axes  fixed  as  in  reference  8;  the  form  of  the  equations  is 
unaffected . 
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k*x 


kl^l  k2^2  k5^5 

*1  +  l2  +  h 


The  right-hand  aide  may  he  expressed  as  x*£,  vhere 

v  _  5k  ^  % 

-  "  l±>  l2>  z3 


defines  a  new  -wave -number  vector • 

Yelocity  at  downstream  station,  general  case,  -  Tlie  velocity  dis¬ 
tribution  whose  curl  in  the  Z^>  £2,  £3  system  is  given  by  equation  (lib) 

and  which  satisfies  continuity  is  found  to  be  expressible  in  the  form 


with 


~  B  *  B  ix-g 
0a  «  Qa  e  “  “ 


0  B  _  1_ 

Hcl  “  V" 


a 


0aA 


-E 


^pAk(3ka 

7  V 
Zp  X 


(12) 


(13) 


where  x  is  the  magnitude  of  the  wave-number  vector  x.  This  result 
is  the  general  solution  for  the  contraction  effect  on  a  single  plane 
wave. 


Equations  (12)  and  (13)  admit  of  a  simple  (but  not  obvious) 
geometrical  interpretation:  traversal’ of  the  stream  contraction  alters 
the  initial  plane  wave,  equation  (8),  so  that  its  wave-number  vector 


kl  k2 

k  =  (k^,  k2,  k^)  is  transformed  into  x  =  — }  ,  —  and  its  amplitude 

- -*r  7 r - X5 — r - w — 7rr .  .  - 

vector  (Q-^ ,  Q2 ,  Q3  )  is  transformed  into  the  projection  of 

Q2A  /z2,  Q^V^)  on  a  plane  normal  to  the  new  wave  -number 
rector  x. 


Yelocity  at  downstream  station,  axis^mmetric  contraction.  -  In 
case  the  stream  contraction  is  axi symmetric*1-  a  considerable  simplifica- 
tion  results.  The  condition  for  axisymmetry  12  *  1 3,  with  use  of  the 

continuity  equation  (2,2),  reduces  equations  (13)  to 


A  contraction  ^uch  that  all  cross  sections  of  the  tunnel  are 
similar,  whence  Z2  =  2 3,  is  termed  axisymmetric;  the  sections  need, not 

be  circular. 
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;  A  ,  2  i  2  ,2 

Ql  +  k2  +  k3 

h  ekx2  +  k22  +  k32 


0  A  QjkjkgCl-e)  ~| 

.  %2  -  *s2  -  *32  J 

3  A  Qikik3(!-e)  ~| 

3  +  ek  2  +  k  2  +  k  2 


> 


J 


where  e= 


(14) 


The  considerably  greater  complexity  of  equation  (13 )  ia  perhaps  obscured 
by  the  purposely  expanded  form  of  equation  (14). 


If  the  initial  wave  normal  k  is  perpendicular  to  the  (longitudinal) 
x-^-axis,  the  component  k-^  vanishes  and  equations  (14)  reduce  further  to 


%B  =  %AAi 
^2B  =  ^2A/J2 
q3b  =  q3a/z2 


(15) 


~  A 

The  same  equations  result  when  Q-j_  may  be  neglected  in  comparison  with 
Q-A  and  that  is,  when  the  amplitude  vector  is  substantially  normal 

to  the  x^-axis.  Equations  (15)  state  that  an  axisymmetric  contraction 


defined  by  Z2 


alters  these  waves  by  a  factor  of 


_1_ 

ii 


in  the  longitudinal  velocity  component  and  a  factor  of  —  in  the  lateral 

2 

velocity  components.  These  equations  apply  only  to  particular  types  of 
waves;  yet  when  the  contraction  effect  is  later  integrated  over  the  random 
aggregation  of  waves  representing  isotropic  turbulence  the  over-all 
results  are  found  not  to  differ  greatly  from  the  simple  factors 

--  and  respectively. 

ll  l2 


The  same  factors  were  obtained  by  Prandtl  (reference  l)  for  other 
special  disturbances:  the  factor  —  from  energy  considerations  for 
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b 


purely  longitudinal  disturbance  velocities,  and  the  factor  -y-  from  momen 

l2 

turn  considerations  for  a  rotating  cylindrical  element  of  fluid  with  axis 
normal  to  the  stream. 


Effect  of  Contraction  on  Spectrum  and  Correlation  Tensors 

Effect  of  contraction  on  correlation  tensor.  -  The  analysis  herein 
leads  first  to  the  changes  in  the  spectrum  tensor  rap^(k)  -*■  rapB  (x ) .  ■ 

Then  the  corresponding  changes  in  the  correlation  tensor  may  he 
obtained  from  the  Fourier  transform  relationship,  equation  (6b): 

GO 

=JjJV(*>e-^(*)  (16a) 

-00 


=, 


RccfP  (jl)  afP  (* ) e  — dT  (i) 


(16b) 


In  succeeding  paragraphs  ra(3B(x)  will  be  determined  in  terms  of  the 
initial  spectrum  tensor  TapA(h)  for  various  cases. 


Spectrum  tensors  at  upstream  and  downstream  stations  in  ter^s  of 
the  Qq.  -  In  an  earlier  discussion  the  Fourier  coefficients  were 

chosen  so  as  to  define  a  field  of  turbulence  confined  to  a  large 
parallelepiped  of  volume  t ,  and  vanishing  everywhere  outside;  for  this 
Lini  8tc^ 

case  <p^00  T  QaQp*  was  to  be  identified  with  the  correlation  tensor 
rap*  I’or  station  A  upstream  of  the  contraction  it  will  be  convenient 

to  specialize  this  parallelepiped  to  a  cube  of  edge  D.  Such  a  cube 
will,  however,  be  distorted  into  a  parallelepiped  of  edges  2-jD,  22D, 

2 3D  by  the  stream  contraction  by  the  time  it  reaches  station  B  down¬ 


stream.  (See  fig.  1(b).)  The  spectrum  tensors  for  stations  A  and 
B,  respectively,  are  therefore 


V®  -^§-QaAte)<vB*© 


B* 


(17) 
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Evaluation  of  spectrum  tensor  at  downstream  station,  general  caseo  - 
The  identifications  made  in  the  last  paragraph  allow  the  spectrum  tensor 
to  he  evaluated  at  station  B  in  terms  of  the  spectrum  tensor  at 
station  A  and  the  parameters  l- Zg,  1 3  defining  the  stream  contrac¬ 
tion  between  stations  A  and  B.  For  a  single  plane  wave  qaA=s  , 

which  is  transformed  by  the  contraction  into  q^  -  equa¬ 

tions  (13)  give 


®  1  n  A 


y  %v. 

t- - J  o  o 


11  *<*[_  r 

In.  the  Fourier  integral  qaE  da  't:o  ^e  interpreted  aa  dq^,  q^®  as 
dlaB^  ^aA  as  QaAdTW>  and  QaB  as  Q^dT  (x) .  Accordingly 


?  -  ^  lb-  E 


,  2  2 

*0  x 


since  =  dT(h)/dT(x).  Thus 


7  2-i  2  7  2 

o  BQ  B*  ll  lZ  l5  0  Aq  A* 


Sh 


A  icykp  QpA  Q^Akyk<x  Q7AQSA 

A2  -  ,2*2  +  2^Zr2x4  ^ 


*7Vx4 


The  corresponding  relation  between  the  post -contraction  and  precontrac¬ 
tion  spectrum  tensors  is,  by  virtue  of  equations  (17), 

(19 ) 


r  b(x)  .  hhh  r  A(k) 

rap  W  lalp  rap  W 


r  ay^  00  k^kp  f  (k )  k^ka  T  (k)  k^kskakp 


,  2  2 
ly  X 


z/x2 


_  Aj  «,  Ll  * 

*7  Z6  x 


where  k  is  related  to  x  by 


kl>  k2>  k3  =  ^lxl>  ^2X2>  J3X3 


Special  case;  axisymmetric  contraction  but  arbitrary  initial 
spectrum.  -  When  the  contraction  is  axi symmetric  =  2^ )  >  th©  equat i on 

of  continuity  in  the  form 
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76 


may  be  used  to  simplify  equation  (19).  The  result  may  be  written 


ap  ^ 


r„*  d)  + 


[ral\  +  rip^aKd-e)  ruS2Vfl(l-*)2 


ek-j_2  +  kg2  +  k32 


ekj2'  +  k22  +  k32 


2  /  2 

where  the  ratio  e  =  Zg  /Z]_  1  f°r  a  large  speed  gain  in  the  contraction 
e  « 1 . 

Special  case:  axisymmetric  contraction,  isotropic  initial  spectrum, 
A  further  simplification  occurs  when  the  turbulence  at  station  A  is 
isotropic.  -  In  that  case,  the  spectrum  tensor  raR®(x)  downstream  of 
the  contraction  can  be  expressed  explicitly  to  within  an  unknown  multi¬ 
plicative  factor  G(k).  This  results  from  the  fact  (reference  3)  that 
r (k )  must  then  be  an  fsotropic  second-order  tensor;  the  isotropic 
property  together  with  the  continuity  condition,  equation  (21),  requires 
Tap  (k)  to  be  of  the  form 

rapAQ£)  =  G(k)(k26^  -  kakp)  (23) 


where 


»  _  f 1  for  a  =  p 

aP  \_0  for  a  /  p 

The  right-hand  side  of  equation  (22)  may  be  evaluated  by  means  of 
equation  (23).  The  diagonal  terms  reduce  to  relatively  simple  forms: 


(*)  = 


(k2  - 

(£k!2  +  k22  +  k32; 


2,  2  /,  2 


r22B,„  .  w  0W  k3 .  ^ H  *-«>  ,  -  *1  m-.> 

Z2  L  ekl  +  k2  +  k32  (ek^2  +  k22  +  k32j 
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-R  T> 

and  r  (x)  is  obtained  from  r22 (*)  by  replacing  k2  by  k3  and 

vice  versa.  The  relation,  equation  (20),  between  k  and  x  applies 
here . 


One -dimensional  longitudinal  spectrums.  -  If  the  form  of  the  initial 
spectrum  tensor  rcxpA(k)is  known,  the  corresponding  one -dimensional 
spectrums  FaA(k-jJ  can  be  calculated,  according  to  the  defining  equa¬ 
tion  (7)  as  applied  at  the  upstream  station  A: 


k2,  k3)dk2dk3 


(26) 


A 


#' 


A  particular  case  of  isotropic  turbulence  is  of  special  interest 
(reference  4):  in  equation  (23)  for  the  function  Gr(k)  is 

taken  to  be  N(k^  +  y2)-3^  ^here  N,7  are  constants.  Then 


3 


and  after  integration 


This  one -dimensional  longitudinal  spectrum  is  of  the  same  form  as  an 
empirical  relation  obtained  in  reference  5  for  that  of  isotropic 
turbulence  in  the  initial  period;  this  agreement  is  the  special  virtue 
of  the  form  assumed  for  G-(k)  • 


The  one -dimensional  lateral  spectrum  functions  corresponding  to 
the  same  G-(k)  are  readily  evaluated;  they  are 


F2A(kx)  =  F3A(k1)  = 


N^k-j2  +  72) 

k2  *  kf 


The  equality  of  the  Fg  and  the  F^  functions  results,  of  course,  from 
the  isotropy  of  the  turbulence. 


The  effect  of  the  stream  contraction  on  these  one-dimensional 
spectrums  is  found  by  employing  the  post-contraction  value  of  F^, 
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tliat  is,  the  -value  appropriate  to  the  downstream  station  B. 

Since  is  a  function  of  the  local  wave-number  vector  x  at 

station  B,  the  equation  corresponding  to  (26)  is 


00 

■«B  - 2  ss 


raa'(Xj/*2,  x  )dx  ax’ 


which,  is  a  function  of  X-j_. 


For  performing  the  integration  and  making  later  comparisons  of 
spectrums  it  is  convenient  to  transform  from  x-j_,  x2,  X3  to  k-]_,  k2, 

k3^  where 

xi  -  ki/h 


x2  =  ^z/l2 

x5  = 


A®  A® 

and  to  define  FaB(k-j_)  =  2]__1  FaB,  such  that  \  Fa®(k1)dk1  =  f  FaBdx^; 


I  r 


fa  k2 
W  h' 


dkgdkg 


The  spectrum  tensor  elements  T  ^  following  an  azi symmetric  contraction 

have  been  evaluated  in  equations  (24)  and  (25).  With  these  values  inserted 
and  Gr(k)  specified  as  before,  the  integrations  of  equation  (25)  are  best 
effected  in  polar  coordinates.  The  results  are  expressed  most  simply  in 
terms  of  a  "normalized "  longitudinal  wave  number  k  as  incorporated 

in  the  two  parameters 

k  2 

- 

a  =  -5—  +  1 

r 

ek-,2  -  kn2 


The  final  result  for  the  one -dimensional  longitudinal  spectrum  following 
an  axisymmetric  contraction  (Z2  =  Z3)  is 
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FlB(kl)  = 


2jtN 


2  2,3 


P*P* 

- 

3  +  4t  +  t2  +  ±_*  +  ( 
2s  V 

2  +  4s  +  2t  +  st  +  -^^logQ 

(a) 

- 

The  corresponding  result  for  the  one -dimensional  lateral  spectrums 
following  an  axi symmetric  contraction  is 


(28) 


(29) 


and  for  e«l  (large  speed  gain)  a  simple  but  very  close  approximation 
is 


F2B(V  =  - 

2 lzZ72 


2  /J2 


i  +  g 

JTT^JTj 


+  0(b,  c  In  e) 


(The  corresponding  approximation  for  F^ip  is  not  simple  enough  to 

warrant  noting.)  The  parameters  Z Zg  =  1 3  and  e  =  1<z‘ [\\  are 

related  to  the  initial  and  final  Mach  numbers  of  the  main  stream  by  the 
equations 


*l‘ 


( mbY  ('■ 5  + 

Vma/  \  5  + 


M, 


'N 


2  %  /  5  +  Mg' 

%  l  5  +  MA2 


Mg' 

3 


h  *  — 


/MaV 
"  W 


5  +  Mb‘ 
5  +  M,  ‘ 


A  /  J 

For  incompressible  flow  (Mg,  MA  -*■  0)  these  reduce  to 

*1  =  VUA 


»2  -  h'l/2 


6  =  Z 


1  J 


(30a) 


(30b) 
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These  post-contraction  spectrums,  equations  (28)  and  (29),  are  com¬ 
pared  with  the  initial  spectrums  in  figures  2  and  3,  respectively;  the 
comparison  is  based  on  an  assumed  initial  stream  Mach  number  of  0,05 
(station  A)  followed  by  an  axi symmetric  contraction  such  that  the  final 
Mach  number  is  2.0  (station  B);  the  corresponding  parameters  are 
-  29.8,  Zg  =  0.382;  e  =  0.00016.  Consider  first  the  longitudinal 

spectrums,  figure  2.  Normalizing  factors  are  incorporated  such  that 
the  areas  under  the  two  curves,  if  replotted  on  a  linear  scale,  would 
be  the  same;  this  normalization  serves  to  differentiate  changes  in 
shape  from  changes  in  amplitude.  The  figure  exhibits  a  rather  striking 
distortion  of  the  spectrum  after  traversing  the  stream  contraction; 
the  peak  spectral  density  is  shifted  from  zero  wave  number  to  kg/y  =  1.4 

along  with  a  general  shift  of  density  to  the  higher  wave  numbers .  Asso¬ 
ciated  with  this  change  in  shape  is  a  reduction  in  amplitude  by  the 

factor  Ug2/uA2,  uA2  and  Ug2  being  the  respective  integrals  of  the 

spectral  density  curves.  These  integrals  are  evaluated  in  a  later 
section. 


The  corresponding  comparison  for  the  lateral  one-dimensional 
spectrums  is  made  in  figure  3.  In  this  case  the  axisymmetric  contrac¬ 
tion  has  made  very  little  distortion  in  the  spectrum.  There  is  again 
a  change  in  magnitude  (this  time  an  increase)  in  the  ratio 

The  changes  in  magnitude  (that  is,  the  changes  in  area  under  the 
spectral  density  curves)  correspond  to  the  changes  ug2/uA2  and 
tb2/va2  in  the  mean  square  components  of  turbulence  and  are,  at  least 

qualitatively,  veil  known.  The  predicted  changes  in  the  shape  of  the 
spectrum  curves  are  apparently  nev„ 

In  the  above  comparisons  both  pre -contraction  and  post -contraction 
spectrums  have  been  expressed  in  terms  of  the  pre-contraction  longi¬ 
tudinal  wave  number  k^,  vhereas  the  local  post-contraction  vave  number 

is  x^  =  Consider,  hovever,  a  representative  longitudinal  vave 

vhich  has  the  form  cos  k-^x  at  station  A  and  cos  at  station  B. 

If  x  and  £  are  identified  vith  the  respective  distances  svept  in 
time  t  by  the  moving  vaves  over  stationary  hot-vire  probes  at 
stations  A  and  B  respectively,  then  k-jx  =  k^Ugt  arwi  x]_£  5=1  xl%^  ” 


—  ZxUAt.  Thus  the  (temporal)  frequency  seen  by  the  hot  vire  in  both 


cases  is  ^1^/2 jc.  The  comparison  based  on  kj_  therefore  constitutes, 
in  effect,  a  comparison  of  the  time  spectrums  that  vould  be  seen  by 
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stationary  hot-wire  probes,  in  contradistinction  to  the  space  spectrums 
discussed  in  the  earlier  parts  of  the  paper. 


Effect  of  Contraction  on  Mean  Square  Yelocity  Components 
for  Isotropic  Turbulence 

The  mean  square  Telocity  components  of  the  turbulent  field  may  be 
identified  as  the  diagonal  terms  of  the  correlation  tensor  Rap(il)  with 
r  set  equal  to  zero.  Thus 

U2  =  B;q(0) 

v2  =  E22(°) 

■W2  =  ^33(0) 


where  u,  v,  w  have  been  written  for  q^,  qg,  qj,  respectively.  The 

evaluation  of  these  means  is  much  less  laborious  than  the  evaluation 
of  the  general  correlation  tensor.  In  particular,  the  evaluation  of 

the  ratio  of  the  means  Ug^/u^,  etc.,  may  be  made  when  the  initial 

turbulence  is  specified  to  be  isotropic  but  no  further  details  of  its 
spectrum,  are  known.  These  ratios  will  be  calculated  in  the  following 
paragraphs . 


Evaluations  of  u2  and  v2  at  upstream  station.  -  According  to 


equation  (6) 


OO 

EnA(0)  .  ^2  .  SJ'J'hiH *><*(*> 


For  isotropic  turbulence  T  ap  has  the  form  specified  in  equation  (23), 


whence 


V 


r\  00 

=jxr< G(k)  (k2  ■  ki2)^) 


where  G-(k)  is  an  arbitrary  function.  It  is  convenient  to  transform, 
to  spherical  polar  coordinates: 
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Because  of  the  unspecified  function  G-(k)  it  is  convenient  to  change 
the  variables  of  integration  from  the  components  of  x  to  the  components 
of  k.  In  other  words,  a  transformation  is  made  from  the  "wave -number 
space"  of  station  B  to  the  ’Hmve-number  space"  of  station  A,  The  trans¬ 
formation  follows  from  the  Cartesian  relations 


dT  (k) 
dT(x) 

together  with  Z2  =  ^3  f°r  an  axisymmetric  contraction,  whence 


=  dk-^dkgdkij 

=  dx-^dxgdx^ 
dk]_  dk2  dkj 
H  l2  l3 


(33) 


Again  the  polar-coordinate  transformation  (equation  (31))  is  made,  with 
the  result 


ub2  =rs 


••  2it  it 

J  k*0(t)  dk  S 


sin 30  d0 


'q  (e  cos^0  +  sin^0) 


The  first  two  integrals  occur  also  in  u^2  (equation  (32)),  and  they 
cancel  in  obtaining  the  ratio  u^2/u^2;  thus 


X  =  _3 

UA2  4^1  v/0 

The  final  result  may  be  written 
2 


sin°0  d6 


(e  cos20  +  sin20) 


u- 


B 


UA 
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T—  + - -T-rr  tanh"1  Yl- 6 

1-6  (1-e)3/2 


(34) 


* 


and  an  asymptotic  expansion  for  small  e  is 
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1  +  —  c  +  (l  +  e)ln  ^  +  0  (£2 


In  e 


Equation  (34)  gives  the  ratio  of  the  mean  square  longitudinal 
velocity  fluctuation  downstream,  of  an  axisymmetric  tunnel  contraction 
to  the  corresponding  mean  square  upstream  of  the  contraction,  when  the 
initial  turbulence  is  isotropic.  The  contraction  is  characterized  by 
an  increase  in  the  stream  speed  in  the  specified  ratio  and  a 

decrease  in  the  lateral  dimensions  in  the  specified  ratio  Z2;  the 
parameters  l±,  Z2,  and  6  =  ll\  are  completely  defined  by  the 
initial  and  final  Mach  numbers  of  the  stream  according  to  equations  (30). 

The  variation  of  */uB2/u/  with  the  speed  ratio  ZB  is  plotted 

in  figure  4  for  two  examples;  in  the  first  the  flow  is  assumed  com¬ 
pressible  with  a  Mach  number  0.05  at  the  start  of  the  contraction;  in 
the  second  the  flow  is  assumed  incompressible  (Ma,  Mg-^0).  The  Mach 

number  scale  at  the  bottom  applies  only  to  the  compressible  case,  the 
l]_  scale  to  both  cases.  The  salient  characteristic  of  the  curve  is 

the  marked  reduction  in  the  longitudinal  component  of  turbulence  with 
increasing  speed  ratio  Z^» 

Compressibility  is  seen  to  have  but  a  secondary  effect,  which  is 
appreciable  only  at  supersonic  speeds.  Note  (equations  (30)  and  (34))  - 
that  with  Z-L  as  the  independent  variable,  the  effect  of  compressibility 

appears  only  in  the  parameter  e .  The  physical  significance  of  e 
follows  from  the  definition  of  Z ^  as  the  speed  ratio  provided  by  the 

contraction  and  Z22  as  the  area  ratio  of  the  contraction  (in  the 
axisymmetric  case  considered),  with  e  =  For  supersonic  final 

speeds  it  is  more  proper  to  speak  of  a  converging-diverging  nozzle  than 
a  contraction,  the  term  "contraction"  having  been  retained  herein 
primarily  for  reasons  of  past  usage. 

The  basis  of  the  compressibility  effect  may  be  summed  up  in  the 
following  way.  The  influence  of  an  axisymmetric  stream  contraction 
arises  from  distortion  of  the  fluid  elements,  as  described  by  the 
parameters  1-^  and  Z2.  (See  fig.  1(b).)  These  parameters  are 

related  by  the  continuity  condition  aZ^g2  =  where  a  is  the 

density  ratio.  Thus  compressibility,  in  allowing  a  to  deviate  from 
unity,  changes  the  relation  between  ZB  and  Z2  somewhat,  and  con¬ 
sequently  modifies  the  contraction  effect. 
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The  graph  of  equation  (34)  in  figure  4  is  primarily  for  illustrative 
purposes;  a  form  more  useful  for  engineering  applications  is  given  in 

figure  5.  The  single  curve  provides  the  variation  of  aJ \ ug2/u^2  with 
"both  and  e;  and  e  may  he  determined  from  the  initial  and 

final  Mach  numbers  by  means  of  the  simple  relation  (30a)* 

p  p 

Evaluation  of  ratio  of  y  at  downstream  station  to  r  at 
upstream  station,  -  The  value  of  v-g^  results  from  an  integration 
involving  the  transformed  spectrum  tensor,  according  to  equation  (16b) 

oo 

e22b(°)  *  VB2  =JJJ r22B(x)dT(x) 

— oo 

Eor  isotropic  initial  turbulence  and  an  axi  symmetric  contraction  the 
transformed  spectrum  tensor  r22B(x)  has  been  evaluated  in  equation  (25). 

Thus 


2kx2k22(i-e) 
ekx2  +  k22  +  k32 


+ 


kx2k22(k2  -  kx2)(l-s)2 
(ekx2  +  k22  +  k32)2 


dT(x) 


Again  it  is  convenient  to  transform  from  x~sPace  to  k-space  (equa¬ 
tions  (33))  and  to  introduce  polar  coordinates  k,  <p,  and  6  (equa¬ 
tions  (3l) )  •  The  integrations  with  respect  to  k  and  <p  are  readily 
disposed  of,  with  the  result 


ain^@  d0  -  2rt(l-e) 


sin^e  cob^O  d 6 
sin20  +  e  cos20 


+ 


«(l-e)2 


ain^e  coa^O  dfl 
aln20  +  e  cos20 


where  H  s 


and  dividing  by 


k4G(k)dk,  as  before.  Upon  carrying  out  the  integration 

o*  Q 

=  —  itH  there  is  obtained  finally 


(35) 


* 
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For  small  e  this  has  the  asymptotic  expansion 


29 


va: 


82- 


^2+e+ie^lneJ  +  0(e2) 


Equation  (35)  gives  the  ratio  of  the  mean  square  lateral  velocity 
fluctuation  downstream  of  an  axisymmetric  tunnel  contraction  to  the 
corresponding  mean  square  upstream  of  the  contraction,  where  the 
turbulence  has  been  assumed  to  be  isotropic.  The  variation  of 


4 


vb2/va2  with  the  speed  ratio  2]_  is  plotted  in  figure  4,  which 


already  contains  the  graph  of  discussed  earlier;  again  the 

two  cases  are  incompressible  flow  and  compressible  flow  with  an  initial 
Mach  number  of  0.05.  For  Jp^l  and  incompressible  flow,  the  lateral 

component  of  turbulence  is  seen  to  increase  steadily  with  2^,  in 

marked  contrast  to  the  decrease  exhibited  by  the  longitudinal  component, 
The  curve  (of  the  lateral  component)  for  compressible  flow  begins  to 
differ  sensibly  from  the  curve  for  incompressible  flow  for  downstream 
Mach  numbers  above  0.3;  above  sonic  speed  compressibility  is  seen  to 
effect  a  complete  reversal  of  the  curve.  The  over-all  effect  of  com¬ 
pressibility  on  the  contraction  effect  is  thus  much  greater  for  the 
lateral  than  for  the  longitudinal  component  of  the  turbulence. 

The  graph  of  equation  (35)  in  figure  4  is  primarily  illustrative; 
a  form  more  useful  for  engineering  applications  is  given  in  figure  6. 

The  single  curve  provides  the  variation  of  aAb^Aa^  with  both 
and  e;  and  e  may  be  determined  from  the  initial  and  final  Mach 

numbers  by  means  of  equations  (30). 


DECAY  CONSIDERATION S 
Criterions  for  Negligibly  Small  Decay 

The  basis  of  the  present  analysis  of  the  contraction  effect  is 
embodied  in  equations  (lla)  relating  the  pre-contract ion  and  post¬ 
contraction  vorticity  distributions.  The  simplicity  of  this  result 
and  its  derivation  arises  from  the  neglect  of  the  turbulent  decay;  by 
decay  is  meant  the  viscous  dissipation  and  all  the  (nonlinear)  inter¬ 
mixing  processes  of  the  eddies  which  together  cause  the  mean  turbulent 
intensity  to  diminish  with  time.  The  postulation  of  an  inviscid  fluid 
eliminated  the  viscous  dissipation,  and  the  limitation  to  very  weak 
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turbulence  eliminated  ttLe  intermixing  processes.  (While  there  can  he 
no  dissipation  in  an  inviscid  fluid,  the  intermixing  processes  ordinarily 
associated  with  decay  will  occur.)  In  order  to  assess  the  influence  of 
these  assumptions  equations  (11a)  will  now  he  derived  in  a  more  general 
fashion  with  the  Navier-Stokes  equations  as  the  starting  point.  For 
simplicity  the  fluid  is  taken  to  he  incompressible,  since  the  major  con¬ 
clusions  are  unaffected  thereby. 

General  formulation  of  changes  in  vorticity.  -  By  rearrangement  and 
cross  differentiation  to  eliminate  the  pressure  term  (reference  8,  p.  578), 
the  Navier-Stokes  equations  can  he  transformed  into 

Dc*L  ^l*  &Li'  &L1'  2 

5T  =  “1^+w2^+“3  +  v?  “l  (36> 

and  two  similar  equations,  where  =00^,002,  003  is  the  vorticity  and 
q*  =  qql,  q.2  %  <13 1  is  the  resultant  velocity,  Now  let  he  the  sum 

of  a  stream  velocity  U,  Y,  W  and  a  turbulent  velocity  field 

£,  ~  I3J  also,  let  curl  U,  Y,  W  =  0,  so  that  ca  is  just  curl 

Then  equation  (36)  becomes  (in  tensor  notation) 


Su  ^1  2 

=  “P  ^  +  “P  ^  +  w  “1 


Contraction  Decay 


(37) 


and  there  are  two  similar  equations.  The  first  set  of  terms  on  the 
right-hand  side  is  identified  as  the  contraction  effect,  the  second 
set  as  the  decay  effect.  First  the  decay  terms  will  he  neglected  in 
an  attempt  to  recover  equations  (lla);  then  the  neglected  decay  terms 
will  he  examined  and  criterions  for  their  neglect  arrived  at. 

Neglect  of  the  decay  terms.  -  Equation  (37)  minus  the  decay  terms 
reads,  in  expanded  form, 


Dco^  d(U+q-]_) 


S(U+qx) 

2  ST  + 


^(U-hi-l) 


In  this  and  the  earlier  equations  —  is  the  fLagrangianT  derivative 
following  the  fluid  motion.  Now  consider  a  line  segment  dx-^,  Sxg, 
following  the  fluid  motion:  its  Lagrangian  derivative  can  he  shown  to  he 
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DSx-j- 
Dt  ' 


S(U+cn) 
1  dXl 


S(u+^i) 

~~^2~ 


+  Sx3 


5(U+ll) 

dx3 


(39) 


and  two  similar  equations.  It  can  be  seen  that  a  solution  of  equa¬ 
tions  (38)  and  (39),  together  with  their  companion  equations,  is  given  by 


®1,  ^2 j  ~  Sxl^  5x2>  8;x3 


for  all  time  t;  this  result  is  well  known.  Now  complete  the  neglect  of 
the  decay  terms  by  omitting  the  terms  in  in  equation  (38)  and  corres 

pondingly  in  equation  (39).  By  this  neglect  the  turbulent  perturbations 
of  the  flow  streamlines  have  been  suppressed:  this  can  be  inferred  from 
the  revised  equation  (39).  If  the  particles  are  at  station  A  at  a  time 
t  =  0  and  reach  station  B  at  time  t  =  t,  there  results 

coqB  Sx-^ 
coxA  5XlA 

Sxq5  Sx-^ 

and  two  similar  equations.  But  - -  is  just  2]_,  - -  is  In, 

5x? 

6X3® 

- -  is  Z3.  Therefore  equations  (lla)  have  been  recovered  for  the 

incompressible  case  (density  ratio  a  =  l) . 


Consideration  of  inertial  decay  terms.  -  In  equations  (37)  to  (39) 
the  decay  terms  not  involving  v  are  the  inertial  or  intermixing  terms. 
These  are  seen  to  be  nonlinear.  The  condition  for  their  neglect  is 
evidently 


GO 


« 


(40) 


and  two  similar  conditions  between  q2  and  Y,  q^  and  W,  respectively, 
In  a  contraction  like  that  of  a  wind  tunnel  the  dominant  velocity 

Su  5y  Sw 

gradients  will  be  and  these  will  be  of  the  same  order  of 

absolute  magnitude.  A  sufficient  condition  to  replace  (40)  is  therefore 


dU 

« 

dx 

(41) 


that  is,  all  of  the  turbulent  velocity  gradients  are  very  much  less 
than  the  axial  gradient  of  the  stream  velocity.  This  is  essentially  the 
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assumption  underlying  the  distortion  equations  (10) ,  which  led  directly 
to  the  vorticity  changes  (lla)  in  Taylor’s  method. 


In  statistical  terms  an  approximate  inference  from  equation  (4l) 
is,  for  isotropic  turbulence, 


But  "by  definition  of  X 


this  may  he  written 


(42) 


The  Microscale1  X  may  he  interpreted  as  a  sort  of  average  eddy  diameter 
weighted  in  favor  of  the  smaller  eddies*  Equation  (42)  may  he  accepted 
as  a  practical  criterion  for  the  neglect  of  the  inertial  decay  terms, 
equivalent  to  one  of  the  two  assumptions  underlying  equation  (10).  The 
other  assumption,  neglect  of  viscosity,  is  considered  next. 


Consideration  of  viscous  decay  term.  -  The  viscous  decay  term  in 
equation" "(37) "  is  the  term  containing  v7  This  term  is  linear  and  so  will 
affect  individual  plane  waves  separately  without  mutual  interference. 

The  magnitude  of  the  term  may  he  estimated  to  a  sufficient  approximation 
hy  considering  a  wave  carried  along  hy  the  contracting  stream  and 
neglecting  (for  this  term  only)  the  distortion  of  the  wave  imposed  hy 
the  contraction.  Thus  a  component  of  the  wave  may  he  written 


„  i(k-x  -  k-,Ut) 

co  i  -  2  ie - 

Then,  if  the  inertial  decay  terms  of  equation  (37)  are  negligible,  the 


equation  reads,  with 

Deo 

Dt 


3u 

du 

^2 

Bx3 

« 

1  5U  2 
_  =  “i  ^ +  vv 

=  “l  +  v  (kl2  +  k22  +  k32)  “1 


I  SU  4.  v2 

=  “1  +  V  y 


* 


% 
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Accordingly,  viscosity  may  "be  neglected  for.  that  portion  of  the  spectrum 
which  satisfies  the  inequality 


vk2  « 


dU 


(43) 


In  the  ’initial  period’  of  decay  if  the  inertial  decay  criterion  (42) 
is  satisfied  the  major  part  of  the  spectrum  will  satisfy  (43). 


Rough  Estimation  of  Mutual  Effects 
of  Decay  and  Contraction 

When  the  decay  effects  are  not  negligible  compared  with  the  con¬ 
traction  effects  (see  criterion  developed  in  the  last  section)  the 
theoretical  basis  of  the  present  theory  of  the  contraction  effect  is 
violated.  Because  negligible  decay  is  more  the  exception  than  the 
rule,  there  is  considerable  incentive  to  attempt  to  apply  the  theory 
outside  the  valid  range  by  means  of  assumptions  concerning  the  simul¬ 
taneous  effects  of  decay  and  contraction. 

Suppose,  now,  the  decay  and  contraction  are  considered  to  occur 
alternatively  in  small  steps,  starting  from  isotropic  turbulence.  Each 
stream  tube  is  considered  to  contract  stepwise:  between  steps  there  is 
decay  without  contraction;  at  each  step  there  is  a  sudden  contraction 
without  decay. _ Let  the  change  in  speed  ratio  per  step  be  d2-j_,  the 

reduction  in  u2  due  to  decay  be  (du2)D,  and  the  reduction  in  u2 
due  to  contraction  be  (du2)c.  Express  the  effect  of  decay  in  the 
absence  of  contraction  in  the  form 

==)  =D(ZX),  (44) 


where  is  a  function  of  the  time  of  travel  (decay  time)  t,  and 

the  effect  of  contraction  in  the  absence  of  decay  in  the  form 


The  corresponding  differential  forms  are 


(45) 
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(du  Jd  D'U-l)  ^ 

'  Biiir a!i 

XT 

(46) 

(du2)C  C'UP  _ 

-T  "  cTuT  ali 

Vi 

(47) 

The  assumption  is  now  made  that  equation  (46)  applies  to  the  decay  effect 
per  step  and  equation  (47)  to  the  contraction  effect  per  step,  the  only 

interaction  being  in  the  common  u  .  The  total  effect  per  step  is  then 


2 

u 


C  1  (2]_) 

c VhT 


V'(i±) 

+  Wh) 


d  l1 


whence  upon  integration  the  over-all  effect  is 

=  =  C(Zj_)  D(Z1)  (48) 

2 

UA 

That  is,  if  the  effect  of  contraction  alone  is  expressed  hy  C(21) 
(equation  (45))  and  the  effect  of  decay  alone  "by  D(2p  (equation  (44)), 

then  the  joint  effect  under  the  assumption  is  expressed  hy  the 
product  C(7p  D(jp. 

Equation  (48)  is  intended  to  provide  a  very  rough  adjustment  of 
the  theoretical  contraction  effect  C(Zi)  to  account  for  decay.  This 


t 


^It  is  known  that  in  the  initial*  period  of  decay  ~  -  u2  dt. 

u2 

Equation  (46)  amounts  to  replacing  the  -u2  on  the  right-hand  side  by 

-  (u2 )decay  only*  some  defense  niade  of  thls  approximation,  con¬ 

sidering  the  progressive  deviation  from  isotropy. 


* 
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adjustment  will  t>©  made  in  the  attempt  to  compare  the  theory  with  experi¬ 
mental  results  in  which  the  decay  effects  are  of  the  same  order  as  the 
contraction  effects . 

Equation  (48)  refers  to  the  longitudinal  velocity  component  u;  an 
equation  of  the  same  form  is  obtained  for  the  lateral  component  v.  For 
both  cases  the  function  D(21)  is  taken  to  be  the  right-hand  side  of 

the  empirical  decay  law  (reference  9) 


for  isotropic  turbulence  in  the  initial  period.  The  decay  time  t(21) 

in  the  formula  is  the  time  required  by  a  particle  of  the  main  stream  to 
pass  through  the  contraction,  the  initial  velocity  being  UA  and  the 
final  velocity  2qUA. 


COMPARISON  WITH  EXPERIMENT 

There  appear  to  have  been  no  experimental  investigations  with  which 
to  compare  the  predicted  changes  imposed  by  a  stream  contraction  on  the 
spectrum  of  the  turbulence,  or  on  the  correlation  tensor  of  the  tur¬ 
bulence.  The  available  experimental  data  seem  to  be  limited  to  measure¬ 
ments  bearing  on  the  changes  in  the  root  mean  square  velocity  components. 
These  data  apply,  moreover,  to  conditions  outside  the  proper  scope  of 
the  present  theory  in  that  large  decay  effects  are  present.  The  experi¬ 
mental  data  are  therefore  compared  with  a  crude  extension  of  the  theory 
in  which  the  decay  is  allowed  for  in  first  approximation.  (See  pre¬ 
ceding  section.) 

The  most  extensive  data  are  those  of  MacPhail,  (reference  10)  which 
in  effect  cover  a  range  of  contraction  ratios  from  2-^  =  1  to  =  9.65 

inasmuch  as  measurements  were  made  at  various  stations  along  the  con¬ 
traction.  Isolated  points  for  particular  contraction  ratios  were 
obtained  from  investigations  made  for  other  purposes  by  Dryden  and 
Schubauer  (reference  ll)  and  by  Hall  (reference  12).  Only  those  points 
were  chosen  for  which  the  initial  turbulence  was  indicated  to  be  approx¬ 
imately  isotropic.  In  the  case  of  reference  11,  data  for  the  case  of 
screens  in  the  settling  chamber  were  excluded  because  the  final  tur¬ 
bulence  level  was  sensibly  indistinguishable  from  the  residual'  noise 
level. 
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In  table  I  are  listed,  for  the  three  experimental  arrangements,  the 
parameters  used  in  the  estimation  of  the  decay  factor  (equations  (48) 
and  (49)).  In  reference  10  the  initial  stream  Telocity  UA  and  integral 

scale  of  turbulence  LA  were  given.  In  reference  11  the  value  of  TJ^ 
was  given,  and  the  value  of  vas  taken  to  be  0.05  feet,  the  only 

scale  mentioned;  it  was  not  clear,  however,  whether  this  value  of  scale 
applied  with  or  without  screens.  In  reference  12  the  value  of  UA  was 

inferred  from  collateral  information  and  is  somewhat  uncertain;  the 
scale  was  estimated  from  the  dimensions  of  the  honeycomb.  In  all 

three  experimental  arrangements  the  initial  relative  levels  of  turbulence 
were  specified.  The  decay  time  t  of  the  turbulence  was  computed  as  the 
time  for  a  particle  to  traverse  the  contraction;  the  value  arrived  at 
for  Hall*s  data  (reference  12)  reflects  the  uncertainty  in  the 
assumed  U^. 

Hoot  mean  square  longitudinal  velocity  components.  -  The  comparison 
of  the  theory,  including  estimated  decay,  with  experiment  for  the  longi¬ 
tudinal  component  of  turbulence  is  given  in  figure  7  •  The  theoretical 
curve,  in  each  instance,  is  the  product  of  a  value  computed  for  con¬ 
traction  alone,  neglecting  decay,  (obtainable  from  fig.  5)  and  a 
value  estimated  for  decay  alone  neglecting  contraction.  (See  eq.ua-' 
tions  (48)  and  (49).)  The  agreement  with  MacPhail’s  data  and  with 
Hall !s  single  point  can  be  considered  good.  The  agreement  with  the 
Dryden-Schubauer  point,  on  the  other  hand,  is  poor;  a  slight  improve¬ 
ment  would  result  on  correction  for  the  spurious  contribution  of  the 
noise  background . 

Hoot  mean  square  lateral  velocity  components.  -  Comparison  of  the 
theory,  again  including  estimated  decay,  with  experiment  for  the 
lateral  component  of  turbulence  is  given  in  figure  8.  There  is  com¬ 
plete  disagreement  with  MacPhail!s  data  and  Hall*s  single  point,  and 
on  the  other  hand,  good  agreement  with  the  Dryden-Schubauer  single 
point.  Thus  there  is  the  curious  result  that  MacPhaills  and  Hallls 
data  agree  well  with  theory  for  the  longitudinal  component  and  disagree 
entirely  for  the  lateral  component,  whereas  the  converse  is  true  for 
the  Dryden-Schubauer  data. 

Discussion.  -  The  uncertainty  both  in  the  manner  of  estimating  the 
decay  effect  and  in  the  data  (table  I)  on  which  the  estimate  was  based 
is  still  far  from  sufficient  to  account  for  the  discrepancies  between 
theory  and  experiment  for  the  lateral  component  of  turbulence.  The 
very  large  amplification  found  by  MacPhail  is  particularly  hard  to 
explain.  On  the  other  hand,  the  experimental  data  of  the  several 
observers  show  considerable  disagreement,  especially  when  differences 
in  decay  are  allowed  for.  This  disagreement  would  tend  to  cast  doubt 
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on  the  validity  of  some  of  the  data;  the  disagreement  may  also  he  in 
part  a  consequence,  predicted  hy  the  theory,  of  possible  differences 
of  the  initial  spectrums  from  each  other  and  from  isotropy, 

concluding  remarks' 

The  original  aim  of  this  paper  was  to  provide  a  quantitative 
explanation  of  the  observed  changes  in  the  root  mean  square  velocity 
components  of  the  turbulence  of  a  wind-tunnel  stream  after  passing 
through  the  tunnel  contraction.  The  simplifying  assumption  of  negligible 
decay  was  made  to  make  the  analysis  tractable,  although  the  decay  and 
contraction  effects  are  ordinarily  comparable.  The  analysis  on  this 
basis  disclosed,  in  addition  to  the  above  integrated  effects,  pronounced 
changes  in  the  spectrum  of  the  turbulence.  The  changes  in  the  shape  of 
the  spectral  density  curves,  as  distinguished  from  over-all  changes  in 
amplitude,  would  appear  to  be  considerably  less  sensitive  to  modifica¬ 
tion  by  decay  than  would  the  mean  square  velocity  components.  For  this 
reason,  and  because  such  spectral  changes  have  not  previously  been  dis¬ 
cussed,  the  emphasis  of  the  present  paper  has  been  placed  most  heavily 
on  these  spectral  effects. 

In  particular,  It  has  been  found  that  the  one -dimensional  longi¬ 
tudinal  spectrum  for  isotropic  turbulence  exhibits  a  rather  interesting 
change  in  shape  downstream  of  the  contraction;  the  center  of  gravity  of 
the  curve  of  spectral  density  versus  longitudinal  wave  number  is  shifted 
substantially  to  higher  wave  numbers,  the  resulting  distortion  moving 
the  peak  of  the  curve  well  to  the  right  of  its  initial  position  above 
the  origin.  The  distortion  is  quite  pronounced  and  would  appear  to  be 
readily  amenable  to  experimental  observation. 

The  restrictive  assumption  of  negligible  decay  largely  defeats  the 
original  aim  of  the  paper.  Nevertheless,  for  practical  reasons  an 
attempt  has  been  made  to  provide  a  crude  extension  to  the  theory  in 
which  decay  is  allowed  for  in  first  approximation.  With  this  approxi¬ 
mation  the  theory  has  been  compared  with  experimental  values  of  the 
contraction  effect  on  the  longitudinal  and  lateral  component  root  mean 
square  velocity  fluctuations.  The  agreement  for  the  longitudinal 
component  is  good,  whereas  there  appears  to  be  almost  complete  disagree¬ 
ment  for  the  lateral  component,  the  experimental  data  themselves  being 
in  conflict.  It  is  perhaps  premature  to  attempt  any  general  conclusion. 
For  the  present,  the  theory  as  augmented  by  the  estimated  decay  effect 
may  be  useful  in  wind -tunnel-design  applications. 

It  is  clear  that  the  tunnel  contraction  effect  on  the  components 
of  turbulent  intensity  cannot  be  represented  by  fixed  fractional  changes 
independent  of  the  character  of  the  initial  turbulence.  Instead  the 
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separate  factors  for  the  longitudinal  and  lateral  components  depend 
markedly  on  the  spectrum  of  the  turbulence.  For  initial  isotropy, 
however,  unique  factors  are  predicted  that,  when  decay  is  neglected, 
are  independent  of  the  details  of  the  spectrum. 

Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  August  30,  1951 
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APPENDIX  A  -  SYMBOLS 

The  following  notation  is  used  in  this  report: 

The  subscripts,  1,  2,  3,  refer  to  a  rectangular  coordinate  system 
with  the  1-axis  alined  with  the  axis  of  the  main  flow  and  directed  down¬ 
stream,  the  2 -axis  directed  horizontally,  and  the  3-axis  vertically,, 
Separate  systems  are  used  with  the  origins  at  stations  A  and  B,  respec¬ 
tively.  (See  fig.  1.)  Vector  and  tensor  notations  are  used  inter¬ 
changeably;  for  example,  k  m  k^  =  (k]_,  k2,  kg) ,  where  a  =  1,  2,  or  3, 

designates  a  vector  with  components  k-j_,  kg,  and  kg. 

C(21)  function  defined  in  equation  (45) 

D(21)  function  defined  in  equation  (44) 


L 


longitudinal  macroscale 
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*3 

M 

K 

Q  = 

1  - 

Eapfe) 

Be 

r 

r  = 


s 


t 

t 

U 

U,V;W  = 

X 


stream  velocity  at  station  B  divided  “by  stream,  velocity  at 
station  A  (see  fig.  l) 

stream,  breadth  at  station  B  divided  by  stream  breadth  at 
station  A  (see  fig.  l) 

stream  height  at  station  B  divided  by  stream  height  at 
station  A  (see  fig.  l) 

Mach  number  of  main  stream 

amplitude  of  special  isotropic  spectrum  tensor  (see  following 
equation  (26)) 


Qa  ~  ($i*  Q3)  disturbance  wave  amplitude  vector 

qa  =  (^i*  q3)  disturbance  velocity  vector 

correlation  tensor  (reference  6) 


real  part  of 

magnitude  of  r  *=  -/r12+r22+r32 

ra  =  (r1<?  r2,  r3)  separation  vector  of 

parameter  in  equation  (28) 
parameter  in  equation  (28) 


two  correlated  points 


time 

main-stream  velocity 

q^,  q2,  q^  disturbance  velocity  components 


length  of  wind  tunnel  contraction  (distance  between  sta¬ 
tions  A  and  B) 


x 


used  occasionally  in  place  of  x^_ 


MCA  TN  2606 


41 


x  = 

7 

e 

6C (07 
e 

X 

X  = 

K 

K  = 

I 

a 

o 

a 

V 

9 

a>  = 

8  = 


xa  =  (x±,  x^}  x^)  position  vector  (station  A) 
spectrum  tensor  (reference  3) 

constant  in  special  isotropic  spectrum  tensor  (see  following 
equation  (26))  (7  ~  l/L) 

contraction  parameter  (=  fl^i  see  fig.  l) 
alternating  tensor  defined  after  equation  (9) 


polar  angle  (equation  (31)) 


magnitude  of 


x12+x22+X32) 


xa  =  (xi>  *3)  transformed  wave  number  vector  (station  B) 
(xa  =  kjla) 

magnitude  of  5  |= 

Ka  =  (?i^  £3)  transformed  position  vector  (station  B) 

(see  equation  (10)  and  fig.  l(*b)) 

summation  over  a  for  a  =  1>  2;  3 


stream  density  at  station  B  divided  "by  stream  density  at 
station  A 

a  volume 

viscosity 

azimuth,  angle  (equation  (3l)) 

coa  =  (co1,  C02;  a)3)  vorticity  vector 

8  =  (8-,  8^,  8*)  vector  amplitude  of  vorticity  wave 

GO  X  O 


Superscripts; 

A  measured  in  vicinity  of  station  A,  upstream  of  contraction 

B  measured  in  vicinity  of  station  B,  downstream  of  contraction 

* 


complex  conjugate 
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Subscripts: 

A  measured  in  vicinity  of  station  A ,  upstream  of  contraction 

B  measured  in  vicinity  of  station  B,  downstream  of  contraction 

take  on  values  1,  2,  or  3  and  designate  tensor  quantities 

1,2,3  specific  values  of  a,  p,  7,  or  8 

A  symbol  with  the  mark  ~  above  it  refers  to  a  single  plane  wave. 

A  bar  over  a  symbol  designates  an  average  (usually  a  spatial  average); 
a  bar  under  a  symbol  designates  a  vector. 


2378 


2378 


MCA  TN  2606 


43 


APPENDIX  B 

EQUIVALENCE  OF  SPACE  AND  TIME  AVERAGES  IN  STATISTICALLY  STEADY, 

HOMOGENEOUS  TURBULENCE 

The  definitions  of  statistical  homogeneity  and  statistical  time 
independence  will  first  "be  made  precise.  Let  F(x,y,z,t)  "be  some 
property  of  a  turbulent  field  that  varies  in  time  and  from  point  to 
point ;  thus  F  may  be  the  pressure,  or  any  of  the  velocity  components, 
or  a  correlation  of  velocity  components  at  two  points  of  fixed  separa¬ 
tion,  (x,y,z)  being  one  of  the  two  points.  If,  for  all  choices  of  the 
property  F,  (a)  the  average  of  F  over  a  time  T  is  independent  of 
(x,y, z),  the  turbulence  is  defined  to  be  statistically  homogeneous;  if 
0>)  the  average  of  F  over  a  volume  V  “► 00  Is  independent  of  t ,  the 
turbulence  is  defined,  in  the  sense  used  herein,  to  be  statistically 
steady  or  time-independent.  The  respective  averages  are  supposed  to  be 
approached  uniformly,  in  the  mathematical  sense,  as  T  or  Y, 
respectively,  approach  infinity.  (A  statistically  steady  or  "stationary” 
condition  is  defined  differently  in  the  theory  of  random  processes.) 

It  will  now  be  proved  that  if  the  turbulence  satisfies  the  two 
conditions  (a)  and  (b),  the  time  and  space  averages  defined  therein  are 
equal.  In  this  proof  no  resort  will  be  made  to  the  "ergodic  hypothesis" 
of  statistical  mechanics,  which  leads  to  the  equt valence  of  the  time 
average  and  the  "ensemble"  average.  The' possibility  of  the  joint 
existence  of  the  conditions  (a)  and  (b)  probably  amounts,  however,  to 
just  as  fundamental  an  assumption. 

The  space  average  will  be  made  over  a  parallelepiped  of  edges  a, 
b,  c  and  the  time  average  over  a  time  T,  and  then  a  limiting  process 
will  be  applied.  The  average  of  F  over  both  space  and  time  is  thus 


F 


s,t 


Lim 

a,b,c,T-*» 


"*T  | 

j 

^b  ' 

/0  ^ 

0  ^ 

lo  V. 

Pa 


F  dx  dy  dz  dt 


(Bl) 


Any  order  of  integration  is  permissible,  since  the  integration 
limits  are  constants.  If  the  time  integration  is  performed  first  the 
expression  may  be  written 


* 

F  +  =  Llm  1 

S,t  a,b,c,T— ^ 

-'c  | 

^b 

j* 

0 

°v 

dx  dy  dz 
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By  virtue  of  the  postulated  uniform,  convergence  of  the  time  and  space 
averages  the  operation  may  he  brought  under  the  integral  sign: 


Lim  I 
a^yC-***  ate 


J 


'T  \  1 

F  dt  |  dx  dy  dz  *  t*?  — 

0  J  J 

(B2) 

dx  dy  dz 


where  F^  is  the  time  average  of  F.  But,  by  condition  (a),  F^  is 
independent  of  x,  y,  and  z.  Therefore 


(B3) 


Alternatively,  the  space  integration  and  limiting  process  may  he 
performed  first: 


fv 

'■'C  | 

■'ll  | 

^a  \ 

—  Lim  i 

/  Lim  _1_ 

a  a  a  \  Lim  1 

Tp  —  J- 

s  ^  t  T-+®  t 

1  a^cr*"  a'bc 

F  dx  dy  dz  dt  -  ^  - 

Vj 

0  \  ^ 

0  ^ 

lo  \_ 

0  / 

(B4) 


where  Eg  is  the  space  average  of  E.  By  condition  (b),  F0  is 
independent  of  t;  therefore, 

*s,t  =  *s  (B5) 

Equation  (B3)  and  (B5)  together  state  that 


(B6) 

or  the  space  average,  the  time  average  and  the  space -time  average  are 
all  equal. 
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TABLE  I  -  DATA  POE  ESTIMATION  OP  DECAY 


*1 

UA 

(ft/sec) 

la 

(ft) 

t 

(sec) 

(ft/sec) 

MacPhail  (reference  10) 

1.20 

3.55 

0.012 

0.19 

0.149 

1.60 

3.55 

.012 

.34 

2.55 

3.55 

.012 

.45 

4.90 

3.55 

.012 

.51 

9.65 

3.55 

.012 

.57 

Dryden-Schubauer 

6.6 

6.86 

.05 

1.31 

.114 

(reference  11 ) 

Hall  (reference  12) 

5.2 

1.54 

.025 

1.22 

.046 

2378 


Figure  5.  -  Variation  of  relative  root  mean  square  longitudinal  velocity  fluctuation  with  both  speed  ratio  12  and  area  ratio 

for  axisyimnetric  contraction  and  isotropic  turbulence  at  J-,  =  1. 
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Figure  7.  -  Comparison  of  predicted  axisymmetric  contraction  effect  with  experiment  for  longi¬ 
tudinal  component  of  turbulence,  with  decay  allowed  for  in  first  approximation.  Initial 
isotropic  turbulence  assumed. 
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